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Abstract —In this paper, we introduce a sophisticated path loss 
model into the stochastic geometry analysis incorporating hoth 
llne-of-sight (LoS) and non-Ilne-of-sight (NLoS) transmissions to 
study their performance impact in small cell networks (SCNs). 
Analytical results are obtained on the coverage prohahility and 
the area spectral efficiency (ASE) assuming hoth a general path 
loss model and a special case of path loss model recommended hy 
the 3rd Generation Partnership Project (3GPP) standards. The 
performance impact of LoS and NLoS transmissions in SCNs 
in terms of the coverage prohahility and the ASE is shown to 
he significant hoth quantitatively and qualitatively, compared 
with previous work that does not differentiate LoS and NLoS 
transmissions. Particularly, our analysis demonstrates that when 
the density of small cells is larger than a threshold, the network 
coverage prohahility will decrease as small cells become denser, 
which in turn makes the ASE suffer from a slow growth or even 
a notable decrease. Lor practical regime of small cell density, the 
performance results derived from our analysis are distinctively 
different from previous results, and shed new insights on the 
design and deployment of future dense/ultra-dense SCNs. 

I. Introduction 

Dense small cell networks (SCNs) are considered as the 
most promising approach to rapidly increase network capacity 
and meet the ever-increasing capacity demands in the 5th gen¬ 
eration (5G) systems IT]. However, up to now, most theoretical 
studies on SCNs only consider simple path loss models that do 
not differentiate Line-of-Sight (LoS) and Non-Line-of-Sight 
(NLoS) transmissions [2-5]. The major conclusion [2-5] is 
that neither the number of small cells nor the number of cell 
tiers changes the coverage probability in interference-limited 
fully-loaded cellular networks. Such conclusion implies that 
the area spectral efficiency (ASE) will monotonically grow 
as small cells go dense. An intriguing question is: Does this 
optimistic conclusion still hold when practical LoS and NLoS 
transmissions are considered in SCNs? 

It is well-known that LoS transmission often occurs when 
the distance between a transmitter and a receiver is small, 
while NLoS is more common in long-distance transmissions as 
well as in office environments and in central business districts. 
For a given network environment, when the distance between 
transmitter and receiver decreases, the probability that an LoS 
path exists between them increases, causing a transition from 
NLoS transmission to LoS transmission. 

To the best of authors’ knowledge, up to now, performance 


analysis considering both LoS and NLoS transmissions are 0 
and Q. In 0 , the authors assumed a multi-slope piece- 
wise path loss function. Such multi-slope piece-wise path loss 
function does not fit well with the NLoS and LoS model 
defined by the 3rd Generation Partnership Project (3GPP) 
standards, in which the path loss function is not a one-to- 
one mapping to the distance 0 . In 17 ], the authors treated the 
event of LoS or NLoS transmission as a probabilistic event for 
a millimetre wave communication scenario. To simplify the 
analysis, the LoS probability function was approximated by 
a moment-matched equivalent step function. The single-piece 
path loss model and the proposed step function for modeling 
the transition from NLoS to LoS transmissions are also not 
compatible with the model recommended by the 3GPP [8, 9]. 

In this paper, we use a general path loss model that features 
piece-wise path loss functions with probabilistic LoS and 
NLoS transmissions. Note that the proposed model is very 
general and includes almost all existing models used to capture 
LoS and NLoS transmissions [6-9] as its special cases. The 
main contributions of the paper are as follows: 

• Analytical results are obtained for the coverage prob¬ 
ability and the ASE using a general path loss model 
incorporating both LoS and NLoS transmissions. 

• Using the above results, closed-form expressions are 
further obtained for the coverage probability and the ASE 
for a special case based on the 3GPP standards. 

• Our theoretical analysis reveals an important finding, 
i.e., the ASE will initially increase with the increase of 
the small cell density, but when the density of small 
cells becomes sufficiently large, the network coverage 
probability will decrease as small cells become denser. 
This in turn makes the ASE suffer from a slow growth 
or even a notable decrease. Thereafter, when the small 
cell density is very large, the ASE will then grow almost 
linearly with the network densification. These results are 
not only quantitatively but also qualitatively different 
from previous study results [2-7]. Thus, our results shed 
new insights on the design and deployment of future 
dense/ultra-dense SCNs in realistic environments. 

The remainder of this paper is structured as follows. 
Section [H] describes the system model. Section [III] presents 


our main analytical results on the coverage probability and 
the ASE, followed by their application in a 3GPP case in 
Section|IV] The derived results are validated using simulations 
in Section |V] Finally, the conclusions are drawn in Section IVTl 


II. System Model 


We consider a DL cellular network in which BSs are 
deployed in a plane according to a homogeneous Poisson point 
process (HPPP) $ of intensity A BSs/km^. UEs are Poisson 
distributed in the considered network with an intensity of 
BSs/km^. A^® is assumed to be sufficiently larger than A so 
that each BS has at least one associated UE in its coverage 
area. The distance between an arbitrary BS and an arbitrary UE 
is denoted by r in km. Considering practical LoS and NLoS 
transmissions, we propose to model the path loss associated 
with distance r as O, shown on the top of the next page. 

In ([TJ, the path loss function (r) is segmented into N 
pieces with the A^-th piece (n (r), n € {1,2,...,W}. For 
each Cn (r), Ch (^) '^he n-th piece of the path loss function 
for LoS transmission, (r) is the n-th piece of the path loss 
function for NLoS transmission and Pr^ (?’) is the n-th piece 
of the LoS probability function. In more detail, 

• Cn (r) is modeled as 


Cn (r) 


CL(r)=ALr-“", for LoS 

(Nl ^ ^ for NLoS 


with and Ajj^,n € { 1 , 2 ,...,W} being the path 
losses at a reference distance r = 1 and aj; and 
a^^,n € { 1 , 2 ,...,W} being the path loss exponents 
for the LoS and the NLoS cases in (n (r), respectively. 
In practice, A{;, aj; and are constants obtained 
from field tests [ 8 , 9]. 

• PrJ; (r) is the n-th piece probability function that a 
transmitter and a receiver separated by a distance r £ 
has an LoS path, which is usually a mono- 
tonically decreasing function of r. For convenience, 
{Pr{; (r)} is further stacked into a piece-wise LoS prob- 


ability function as 

'PrL(r), 

when 0 < r < di 


Pr^ (r), 

when di < r < d 2 

Pr^ (r) = < 

■ (3) 



when r > dx-i 


Our model is consistent with the ones adopted in the 
3GPP [ 8 , 9]. Note that the considered path loss model in ([T]i 
will degenerate to that adopted in lb) and Q when Pr{; (r) = 
0, Vn € {1, 2,... ,N} and N = 1, respectively. 

As a common practice in the field [2-6], each UE is assumed 
to be associated with the nearest BS to the UE, and the multi- 
path fading between an arbitrary BS and an arbitrary UE is 
modeled as independently identical distributed (i.i.d.) Rayleigh 
fading, i.e., the channel gain is denoted by h and is modeled as 
an i.i.d. exponential random variable (RV). The transmit power 
of each BS and the additive white Gaussian noise (AWGN) 
power at each UE are denoted by P and Nq, respectively. 


III. Analysis Based on General Path Loss Model 


Using the properties of the HPPP, we study the performance 
of SCNs by considering the performance of a typical UE 
located at the origin o. We first investigate the coverage 
probability and thereafter the ASE. 

The coverage probability is defined as the probability that 
the signal to interference plus noise ratio (SINR) of the typical 
UE, denoted by SINR, is above a threshold 7 : 

p“''(A, 7 ) = Pr[SINR> 7 ], (4) 


where the SINR is computed by 
SINR = 


PC (r) h 


(5) 


( 6 ) 


Ir + Nq ’ 

where Ir is the cumulative interference given by 

Ir — ^ ^ PPi9ii 

where bo is the BS associated with the typical UE and located 
at distance r from the typical UE, and /3i and gi are the path 
loss and the multi-path fading channel gain associated with 
the i-th interfering BS, respectively. 

According to 0 and m, the area spectral efficiency (ASE) 
in bps/Hz/km^ for a given A can be computed by 


^ (A,7o) = ^/ log 2 (1 +t)/x (A,a:)da;, (7) 

J 70 

where 70 is the minimum working SINR for the considered 
SCN, and fx (A, x) is the probability density function (PDF) 
of SINR observed at the typical UE at a particular value of A. 

Since (A, 7 ) can be defined as the complementary 
cumulative distribution function (CCDF) of SINR, fx (A, x) 
can be expressed as 

5(1-p“''(A,t)) 


fx (A, x) = 


dx 


( 8 ) 


Given the definition of the coverage probability and the ASE 
presented in (|4]i and (|7J respectively, and using the path loss 
model of O, we present our main result on p™'' (A, 7 ) in 
Theorem [T] shown on the next page. Plugging p‘^°'' (A, 7 ) from 
0 of Theorem [T] into we can get the ASE using (|7]i. 

As can be seen from Theorem [1] the coverage probability, 
p‘^°''(A, 7 ), is a function of the piece-wise path loss func¬ 
tion {Cn(T)} and the piece-wise LoS probability function 
{PrJ; (r)}. We will investigate their impacts in the sequel. 


IV. Study of a 3GPP Special Case 

As a special case of Theorem [T] we consider the path loss 
function, ^ (r), 

w/ probability Pr^ (r) 

probability (l — Pr^ (r)) 

together with the linear LoS probability function, which is 

' 1-37’ 0<r<di 

0 , r > di 


Cir) = 


Pr^ (r) = 


(14) 


(15) 


both respectively recommended in the 3GPP [ 8 , 9]. 

Considering the general path loss model presented in ([1]), 
the path loss model presented in ([T4l i and dTSl) can be deemed 






c (r) = < 


C 2 (r) = 


Civ (r) = 


C\{r), 

cr(o, 

with probability Pr^f (r) 

with probability (l — Pr^f (r)) ’ 

when 0 < r < di 

C2"(0, 

CT{r), 

with probability Pr 2 (r) 

with probability (l — Pr^ (r)) ’ 

when di < r < d 2 

Cat (r) , 

with probability Pi"*^ (r) 
with probability (l — Pi"^ (r)) 

, when r > djv-i 


( 1 ) 


Theorem 1. Considering the path loss model o/([7J, p™*" (A, 7 ) can be derived as 

p"''(A,7) = 


(9) 


where Til = > 7] fR,n (f) dr, Pr > 7] [r] dr, and do and dN are respectively defined , 

0 and 00 . Moreover, (r) and (^) t^t-e represented as 


arid 


/fl,n {r) = Pfn {r) X exp (-Trr^A) x 27rrA, (d„_i < r < dn) , 

fii,n (r) = (1 - Pr^ (r)) x exp (-Trr^A) x 2-7rrA, (d^-i < r < dn) . 


Pr 


and 


Pr 


[ 7r + iVo ^ 'J 

are respectively computed by 

\PCn{r)h 

> 7 

= exp 

f 7(Vo ^ 


Ir + A^O 

V P(,ii{r)) 

PC{r)h 

> 7 

= exp ^ 

yNo 


Ir + A^O 

Pa^r), 


7 


P^i.{r)J ’ 

7 

PCn"- (r) 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


where T£in (s) is the Laplace transform of RV T evaluated at s. 
Proof: See Appendix A. 


According to Theorem [T] p“'' (A, 7 ) can be obtained as 

2 

(A, 7) = ^(T^ + rr). (16) 


In the following subsections, we investigate T^, T 2 , 

■NL 


and respectively. 

A. The Computation of for 3GPP Case 1 


From Theorem [T] T]f can be derived as 


rdi 

^1 = / 

Jo 


Tf = exp I ^ ] /]i,i (r) dr, (17) 


where according to Theorem [T] and (fTSl l. (r) becomes 

fki w = 

n — ^ j X exp (— Trr^A) x 27rrA, (0 < r < di). (18) 

Furthermore, to compute I -p^ ) in the range of 0 < 
r < di, we propose Lemma |2] 


calculated by 

(. 

1 Ta^ 


as a special case of O with the following substitution: N = 2, Lemma 2. jSf/ ( ^ 1 in the range of 0 < r < di can be 
Ci- (r) = It) = ALr-“ , Cr (^) = (^) = \ pa I 

Prlj; (r) = 1 — and Pr^ (r) = 0. For clarity, this 3GPP 
special case is referred to as 3GPP Case 1 in the sequel. 


1 2 I_ 1 - 


exp ^-27rA 1, (^ 7 r“^^ ,di^ 

-pi 1, (7r“‘') 

^a^, 2 , (^ 7 r“'') ,di^ 


X exp 


V d,o 




X exp (pi 


(l^r° 
V AL ^ 


-1 


, di 


-Pi 


exp ^-277Ap 2 1, ^ 


where 


(0<r<di), (19) 






























Pi {a,P,t,d) = 
■c?(/3+l) 


P 2 {a,P,t,d) = 
t{a — P — 1) 


2 F 1 


2 F 1 




a 


( 20 ) 


’ a ’ a ’ td° 


, ( 21 ) 


where 2 F 1 [•,•;•;•] is the hyper-geometric function m- 
Proof: See Appendix B. 

B. The Computation of for 3GPP Case 1 
From Theorem [1] can be derived as 


rr=r;xp(-2;S^i^, 


TA^ j (^) ( 22 ) 

where according to Theorem [T] and (fTSl l. f^i (r) becomes 
(’") = ^ X 27rrA, (0 < r < di). (23) 

«i 


) 


^NL 

Furthermore, to compute ) in the range of 0 < 

< c?i, we propose Lemma [3] 


,L 1 


^NL 


-pi I a^,l, [ 


X exp 



Th = 




Furthermore, to compute I in the range of r > 


Lemma 3. | range of 0 < r < di can be 

Iculi 

exp —27rA I pi I a^, 1 , 


-PI ^Q^2, I 

X exp (pi 2, , di^ 

-pi 2, (7?'“'''') ’ ^) ) ) 

X exp ^-27rAp2 1, ( 7 ’'“ ) , ^ , 

(0 < r < di), (24) 

where pi (a, /3, t, d) and p 2 (ct, /3, f, d) are defined in ( 1201 ) and 
( 1271 ). respectively. 

Proof: See Appendix C. ■ 

C. The Computation ofT^for 3GPP Case 1 
From Theorem [1] can be derived as 


= 0 . 


(25) 


Note that the reason why Tf = 0 \n (I25l l is because 
according to Theorem [1] and (fTST l. we have 

/fl ,2 (^’) = Ox exp (—Trr^A) x 27rrA 

= 0, (r > di). (26) 

D. The Computation of for 3GPP Case 1 
From Theorem [T] can be derived as 

=/Jexp (0 dr^ ill) 

where according to Theorem [1] and (fTST l. f ^2 (^) becomes 
/h ^2 (^) = lx exp (—Trr^A) x 27rrA 

= exp (—7rAr^) x lirrX, {r > di). (28) 

ijNL ' 

Ir \ P.4NL 

di, we propose Lemma H] 

/ cNL \ 

Lemma 4. I ) in the range of r > di can be 

calculated by 

= 

exp ^- 2 pAp 2 1 , ( 7 ^"'""^) . V > <ii) , (29) 

where p 2 (ex, /3, t, d) is defined in ( 1271 ). 

Proof: See Appendix D. ■ 

E. The Results of p™'’ (A, 7 ) and (A, 70 ) 

To sum up, (A, 7 ) for 3GPP Case 1 can be written as 
p“''(A, 7 ) = T\ + Tf^ + T^^, (30) 

where T^, and are computed from closed-form 

expressions using (fTTl i. (l 22 l i and (l27l i. respectively. 

Plugging (A, 7 ) obtained from (l30l) into (0, we can get 
the result of (A, 70 ) from d?) for 3GPP Case 1. 

V. Smulation and Discussion 

In this Section, we use simulations to further study the 
performance of SCNs and establish the accuracy of our 
analysis for 3GPP Case 1 studied in Section HV] According 
to IS] and ||9l, we use the following parameters: di = 0.3km, 
= 2.09, = 3.75, A^ = lO-^ ^, A^^ = lO-^-^s, 

P = 24 dBm, A^o = -95 dBm. 

A. Validation and Discussion of p™'’ (A, 7 ) 

The results of p‘^°'' (A, 7 ) with 7=1 and 7 = 10 are plotted 
in Fig.[T] As can be observed from both figures, our analytical 
results perfectly match the simulation results. Since the results 
of A'^^^(A, 7 o) are computed based on p^°'^ we will 

only use analytical results for p“'' (A, 7 ) in our discussion 
hereafter. For comparison, we have also included analytical 
results assuming a simplistic path loss model that does not 
differentiate LoS and NLoS transmissions 0 . Note that in 0 , 
only one path loss exponent is defined and denoted by a. Here, 
a is set to to show the analytical results from 0. 

































BS density A. [BSs/km^] 

Fig. 1. (A, 7 ) vs A for 3GPP Case 1. 

From Fig. [T] we can observe that the coverage probability 
performance given by the stochastic geometry analysis in Q 
first increases with the BS density because more BSs provide 
better coverage in noise-limited networks. Then, when A is 
large enough, the coverage probability becomes independent 
of A since the network is pushed into the interference-limited 
region, e.g., A > 10^ BSs/km^ for the analysis in Ha. This 
observation is consistent with the conclusion in El, which 
shows that for a sufficiently large A, the coverage probability 
becomes almost a constant with the increase of the small cell 
density. The intuition behind the observation is that with the 
simplistic assumption on the path loss model, the increase in 
interference power is counterbalanced by the increase in signal 
power in a interference-limited network, and thus the coverage 
probability remains the same as A increases. 

In Fig. [T] the coverage probability performance of the 
proposed stochastic geometry analysis incorporating both LoS 
and NLoS transmissions exhibits a signihcant deviation from 
that of the analysis from El, because when the distance 
r decreases, or equivalently when the small cell density A 
increases, LoS transmission occurs with an increasingly higher 
probability than NLoS transmission. Specihcally, when the 
SCN is sparse and thus noise-limited, e.g., A < lOBSs/km^, 
the coverage probability given by the proposed analysis grows 
as A increases for the same reason as explained in the above 
paragraph, i.e., deploying more small cells is beneficial for 
removing coverage holes. Then, when the network is dense 
enough and all coverage holes are removed, the coverage 
probability given by the proposed analysis decreases as A 
increases, due to the transition of a large number of interfer¬ 
ence paths from NLoS to LoS. It is important to note that the 
coverage probability performance of the proposed analysis for 
3GPP Case 1 peaks at a certain value Aq. When A increases 
above Aq, interfering BSs may be very close to the typical 
UE and hence their signals may reach the UE via strong LoS 
paths. Such critical point of Aq can be readily obtained by 
setting the partial derivative of (A, 7 ) with regard to A 
to zero, i.e., Aq = arg | ~ solution to this 

equation can be numerically found using a standard bisection 
searching IfTTI . In Eig. [ 1 ] the numerical results of Aq are 



Fig. 2. (A, 70 ) vs A for 3GPP Case 1. 

15.85 BSs/km^ and 10.21 BSs/km^ for 7 = 1 and 7 = 10, 
respectively. 

B. Discussion of the Analytical Results of 

The results of 2 ) and (A, 10) are plotted in 

Fig. El comparing the proposed stochastic geometry analysis 
with the conventional stochastic geometry analysis in El- 
As can be seen from Eig. El the analysis from |EI indicates 
that when the SCN is dense enough, e.g., A > 10^ BSs/km^, 
the ASE performance increases almost linearly with A, which 
is logically correct from the conclusion that (A, 7 ) is 
invariable with respect to A for a given 7 in dense SCNs El- 
In contrast, the proposed stochastic geometry analysis re¬ 
veals a more complicated trend for the ASE performance. 
Specifically, when the SCN is relatively sparse, e.g., A < Aq, 
the ASE quickly increases with A because the network is 
generally noise-limited, thus adding more small cells im¬ 
mensely benehts the ASE. When the SCN is extremely dense, 
e-g-, A > lO^BSs/km^, the ASE exhibits a nearly linear 
trajectory with regard to A because both the signal power 
and the interference power are now LoS dominated and thus 
statistically stable. Note that the pace is lower than in El- As 
for the practical range of A, i.e., A £ [Aq, 10^] BSs/km^, the 
ASE first exhibits a slowing-down in the rate of growth (when 
7 = 1 ) or even a notable decrease in its absolute value (when 
7 = 10). This is attributed to the fast decrease of the coverage 
probability at around A £ [Aq, Ai] BSs/km^ as shown in Eig.[T] 
where Ai is another threshold larger than Aq. When A > Ai, 
the ASE will pick up the growth rate as the decrease of the 
coverage probability becomes very gentle. In Eig. El the value 
of Ai seems to be around lO^BSs/km^. 

Our new hnding indicates the signihcant impact of the path 
loss model incorporating both NLoS and LoS transmissions. 
As a conhrmation, noting that in Eig. [T1 we can observe that 
increasing 7 from 1 to 10 will greatly accelerate the decrease 
of the coverage probability at around A £ [lO, 10^] BSs/km^, 
which in turn causes the notable decrease of the ASE at that 
range of A when 7 = 10 in Eig. El 

With the dehned thresholds Aq and Ai, SCNs can be roughly 
classihed into 3 categories, i.e., the sparse SCN (0 < A < Aq), 






















the dense SCN (Aq < A < Ai) and the very dense SCN (A > 
Ai). The ASEs for both the sparse SCN and the very dense 
SCN grow almost linearly with the increase of A, while the 
ASE of the dense SCN shows a slow growth or even a notable 
decrease with the increase of A. Erom Eig.|2] we can get a new 
look at the ultra-dense SCN, which has been identified as one 
of the key enabling technologies of the 5G networks fT]. Up to 
now, there is no consensus in both industry and academia on 
that at what density a SCN can be categorized as an ultra-dense 
SCN. According to our study, for 3GPP Case 1, we propose 
that the 5G systems should target the third category of SCNs 
as ultra-dense SCNs, i.e., the SCNs with A > Ai, because 
the associated ASE will grow almost linearly as A increases. 
Numerically speaking, Ai is around 10^ BSs/km^ from Eig. |2] 
It is particularly important to note that the second category 
of SCNs (Ao < A < Ai) is better avoided in practical SCN 
deployments due to its cost-inefficiency shown in Eig. |2] 


VI. Conclusion 

In this paper, we have shown that a sophisticated path loss 
model incorporating both LoS and NLoS transmissions has 
a significant impact on the performance of SCNs, measured 
by the two metrics of the coverage probability and the ASE. 
Such impact is not only quantitative but also qualitative. 
Specifically, our theoretical analysis have concluded that the 
ASE will initially increase with the increase of the small 
cell density, but when the density of small cells is larger 
than a threshold Aq, the network coverage probability will 
decrease, which in turn makes the ASE suffer from a slow 
growth or even a notable decrease as the small cell density 
increases. Eurthermore, the ASE will grow almost linearly as 
the small cell density increases above another larger threshold 
Ai. According to our study, for 3GPP cases, we propose that 
the 5G systems should target the SCNs with A > Ai as ultra- 
dense SCNs. Numerically speaking, Ai appears to be around 
several lO^BSs/km^. 

The intuition behind our conclusion is that when the density 
of small cells is larger than a threshold, the interference power 
will increase faster than the signal power due to the transition 
of a large number of interference paths from NLoS to LoS, 
and thus the small cell density matters! 

As our future work, we will incorporate more sophisticated 
UE association strategies and more practical multi-path fading 
model into the analysis of SCNs because the multi-path fading 
model is also affected by LoS and NLoS transmissions. 


as 


Appendix A; Proof of Theorem[T] 

Erom dUi and (|5]l, we can derive (A, 7 ) straightforwardly 

p“''(A, 7 ) = f Pr[SINR> 7 |r]/i^(r)dr 

J r>0 

PC (r) h 


I 


Pr 




r>0 

N 


■ Nn 


> 7 


Ir (r) dr 


n—1 


where Tjr and T, 


NL 


^ n 
d. 




d-n- 


= L 

Pr 


dn 


Pr 


are piece-wise functions defined as 


PCCir)h 


dn — 

Ir+No ^ ' 


> 7 


fR,n (r) dr 


and T, 


-NL 


Ir + No 

f^n (^) respectively. Besides, do 
and c?Ar are respectively defined as 0 and 00 . Moreover, 
fh. n (’’) f^n (^) '^he piece-wise PDEs of the event 
that the UE is associated with the nearest BS with an LoS 
path at distance r and the event that the UE is associated with 
the nearest BS with an NLoS path at distance r, respectively. 

Regarding „ (r), we define two events in the following, 
whose joint event is equivalent to the event that the UE is 
associated with a BS with an LoS path at distance r. 

• Event B: the nearest BS is located at distance r 

• Event D^: the BS is one with an LoS path 
According to the cumulative density function (CDE) of 
Event B with regard to r is given by 

PR,n (r) = 1 - exp (-Trr^A) , (d„_i < r < d„). (32) 

Hence, taking the derivative of (r) with regard to r, 
yields the PDE of Event B as 

/fi.n (r) = exp (-Trr^A) x 27rrA, (d„_i < r < d„). (33) 

The PDE (r) should be further thinned by the probability 
of Event on condition of r, which is Pr^ (r), so that we 
can get the PDE of the joint event of B and as 

fR,n (r) = Pr^ (r) X /|„ (r). (34) 

Regarding (r), we also define two events in the fol¬ 
lowing, whose joint event is equivalent to the event that the 
UE is associated with a BS with an NLoS path at distance r. 

• Event B: the nearest BS is located at distance r 
« Event the BS is one with an NLoS path 

Similar to (I 34 I 1 . the PDE fi,nir) should be further thinned 
by the probability of Event on condition of r, which is 
(1 — Pr^ (r)), so that we can get the PDE of the joint event 
of B and as 

(1 — Pr^ (r))xexp (—Trr^A) x27rrA, (c?„-i< r< d„). (35) 


As for the calculation of Pr 

'PCkir)h 


PQ{r)h 
Ir + No P ’ 


Pr 


■ ZVn 


> 7 


= ]E 


[Ir 


= E 


[Ir 


Pr 


Fh 


h > 


in (EB, we have 

7 (Ir + No) 


PCk ir) 

7 {Ir + Nq 

PCkir) 


(36) 


where E[x] {•} denotes the expectation operation by taking 
the expectation over the variable X and Fh [h) denotes the 
complementary cumulative density function (CCDE) of RV h. 
Since we assume h to be an exponential RV, we have Fh (h) = 
exp {—h) and thus (l36l) can be further derived as 
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where .2/ (s) is the Laplace transform of evaluated at s. 
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Our proof of Theorem [D is completed by plugging 0^ . 
ED, dS, and ED into ED- 

Appendix B; Proof of Lemma[2] 

Based on Theorem [D it is straightforward to derive (s) 

in the range of 0 < r < di as 

(s) = E[/^] {exp(-s7r.)|0 < r < di} 
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0<r<dij, (39) 

where (a) in EH* is obtained from 0. 

Since 0 < r < di, E[g] {exp (—sP/3 (u) g)} in (1^ should 
consider interference from both LoS and NLoS paths. Thus, 
(s) can be further derived as 


■^ir (s) = exp ( —27rA 
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+ (sPANL)-iu“"''^ ' 


Plugging s = 7^417 into (l40l i. and considering the definition 
of Pi {a,l3,t,d) and p 2 {a,(3,t,d) in (l20l) and (|2D, we can 


obtain pW J shown in (fT^ . which concludes our proof. 

Appendix C: Proof of Lemma[3] 

Following the same approach in Appendix B, it is ready to 
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derive jSf/,, I 1 in the range of 0 < r < as 
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Similarly, ^ paott 
be calculated by 


in the range of xi < r < di can 
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Our proof is thus completed by plugging (l20l i and (1211) into 
(ED and (ED. 


Appendix D; Proof of Lemma E] 

Following the same approach in Appendix B, it is ready to 

/ £,NL\ 

derive 2?/,, I 1 in the range of r > di as 

^ =exp f- 27 rA [ - - -j- du\ 

= exp ^-27rAp2 ) ,r^^,(43) 

where p 2 (a, /3, t, d) is defined in (l2TT i. 

Our proof is thus completed with (ED . 
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